Head-on collision of compact objects in general relativity: 
Comparison of post-Newtonian and perturbation approaches 
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The gravitational-wave energy flux produced during the head-on infail and collision of two com- 
pact objects is calculated using two approaches: (i) a post-Newtonian method, carried to second 
post-Newtonian order beyond the quadrupole formula, valid for systems of arbitrary masses; and (ii) 
a black-hole perturbation method, valid for a test-body falling radially toward a black hole. In the 
test-body case, the methods are compared. The post-Newtonian method is shown to converge to the 
"exact" perturbation result more slowly than expected a priori. A surprisingly good approximation 
to the energy radiated during the infall phase, as calculated by perturbation theory, is found to be 
given by a Newtonian, or quadrupole, approximation combined with the exact test-body equations 
of motion in the Schwarzschild spacetime. 



PACS numbers: 04.25.Nx, 04.30.-w, 97.60. Jd, 97.60. Lf 



I. INTRODUCTION AND SUMMARY 



The head-on collision of two bodies in general relativity has long been a test-bed for various approaches to solving 
Einstein's equations in dynamical situations, despite its vanishingly small astrophysical probability. Because of the 
■ axial symmetry of the situation, the problem simplifies dramatically, yet can still retain features of potential astro- 
physical importance, such as strong curvature, fast motions, and the emission of gravitational radiation to infinity. 
The first attempt to solve this problem in a fully general relativistic manner was the 1971 calculation of perturbations 
of the Schwarzschild black- hole spacetime caused by a test particle in radial infall from infinity |l]] , using the recently 
developed Zerilli equation for black-hole perturbations || . The main result of this calculation was a formula for the 
total gravitational- wave energy emitted, AE = 0.0104 fic 2 (fi/m), where fj, is the mass of the test body, and m the 
mass of the black hole (fj, -C to). Subsequently, these perturbation methods were extended to handle more general 
>» ■ test-body orbits around black holes (^-[|] ■ 

The head-on collision of two black holes has also been the starting point in attempts to solve general relativistic, 
I . dynamical problems by purely numerical techniques. The pioneering computation of the head-on-collision of two 
^Jy equal-mass black holes in the middle 1970s by Smarr and Eppley || stood for almost 15 years as the state of the 
art, until recent technical advances together with the use of supercomputers resulted in a substantial improvement in 
accuracy and reliability |tJ . Generalizing these computations to the case of black holes in circular orbits is sufficiently 
difficult that it is viewed as a "grand challenge." 

Another method for studying the interaction of compact objects and the emission of gravitational radiation is the 
post-Newtonian approximation. Although this method can be used for bodies of arbitrary masses in arbitrary orbits, 
it breaks down when the gravitational fields become strong and the velocities become comparable to the speed of 
light. 

The ultimate goal toward which these calculations are pointed is the solution of the inspiral and coalescence of 
two compact objects (black holes or neutron stars), and the calculation to high accuracy of the emitted gravitational 
waveform. Such inspiralling systems are believed to be the most promising sources for detection of gravitational 
waves by kilometer-scale laser interferometric gravitational observatories (LIGO ||] in the US, VIRGO Jj| in Europe). 
Furthermore, it has been shown that very accurate gravitational-waveform templates will be needed in order to infer 
characteristics of the sources from the data, using matched filtering techniques |l(J. It is also widely believed that 
all three theoretical techniques will play a role in this effort: post-Newtonian techniques for describing much of the 
inspiral phase (until the system becomes too relativistic), numerical techniques for the late inspiral and coalescence, 
and test-body perturbation techniques for the special case of a small mass inspiralling into a very massive black 
hole. Consequently, it is important to explore carefully the regions of validity and the accuracy of these techniques, 
especially in regimes where they overlap. 

In this paper we study the gravitational radiation emitted during the head-on collision of two compact objects. 
In particular, we compute the gravitational-wave luminosity as a function of their relative separation. We use two 
techniques: (i) the post-Newtonian approximation, carried to second order [order (v/c) 4 ~ (Gm/rc 2 ) 2 ] beyond the 
quadrupole approximation; and (ii) the test-body perturbation technique. Specializing the post-Newtonian method 
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to the case of a test-body falling from infinity toward a massive body, we compare the two approaches. Our main 
findings are these: 

(1) The post-Newtonian approximation begins to fail when the bodies are separated by about lOGm/c 2 (in harmonic 
coordinates), where m is the total mass of the system. By failure, we mean that post- Newtonian corrections to the 
energy flux become comparable to the leading-order, quadrupole term. 

(2) The post-Newtonian approximation converges very slowly: at a separation of 100Gm/c 2 , the first and second 
post-Newtonian corrections to the luminosity are respectively 5 and 40 times their expected sizes of Gm/rc 2 ~ 10 -2 , 
and (Gm/rc 2 ) 2 ~ 10~ 4 , relative to the leading-order term. 

(3) The test-body approach shows that for separations larger than AGm/c 2 (in Schwarzschild coordinates), the 
radiation emitted is dominated by bremsstrahlung radiation generated by the infalling particle. Thereafter, the radi- 
ation is dominated by black-hole quasi-normal ringing. The bremsstrahlung radiation contributes only approximately 
3 percent of the total energy emitted. 

(4) At separations greater than about 10Gm/c 2 , the perturbation results can actually be well approximated by a 
hybrid "quadrupole" formula for the luminosity, in which the exact Schwarzschild test-body equations of motion are 
used to evaluate the three time derivatives of the source quadrupole moment. The post-Newtonian values differ from 
the perturbation results by 20% at a separation of A5Gm/c 2 , 10% at 75Gm/c 2 , and 1% at 400Gm/c 2 , reflecting a 
very slow convergence. The slow convergence of the post-Newtonian approximation has also been noted in the case 
of circular test-body orbits O-H. 

The organization of this paper is as follows: Section II discusses the head-on collision of bodies of arbitrary masses 
to second post-Newtonian (2PN) order beyond the quadrupole approximation, including the effects of gravitational- 
wave tails. Section III discusses the perturbation approach, and presents the energy flux for test-body radial infall 
as a function of retarded time, showing both the bremsstrahlung radiation and the dominant quasi-normal ringing 
radiation. In section IV we compare the two approaches; to do so requires a careful mapping between the post- 
Newtonian retarded source variables, expressed in harmonic coordinates, and the retarded Schwarzschild coordinates 
of the perturbation approach. Section V discusses the accuracy of the post-Newtonian approach, and derives the 
hybrid "quadrupole" formula for the luminosity. In an Appendix we discuss a remarkable series of cancellations that 
occur in the post-Newtonian method for the radial infall from infinity. 



II. HEAD-ON COLLISION OF COMPACT OBJECTS TO SECOND POST-NEWTONIAN ORDER 
A. Gravitational wave generation to second post-Newtonian order 

This first subsection summarizes from the literature the relevant equations for gravitational wave generation by a 
slowly moving system of two arbitrary, non spinning, masses Jl5| . 

We describe the motion of the two bodies in terms of their relative distance. In a harmonic coordinate system whose 
origin is at the center of mass, the relative position vector is x = X\ — X2, where X\ and X2 are the positions of each 
body. If the masses arc denoted m\ and 7712, we define the total mass m = mi +TO2, the reduced mass fi = m\ m^/rn, 
■q = [i/m, and 5m = m\ — m^. Also, r = \x\ and v = x, where a dot denotes differentiation with respect to time. 

The luminosity (or energy flux) associated with the gravitational waves can be expressed in terms of the mass and 
current multipole moments of the radiation field |l6[ ] . We need only a few terms for an expression accurate through 
second post-Newtonian (2PN) order beyond the leading-order, quadrupole-formula expression: 

E < r > = \ ( (311 'f + m ( m * m Y + tea C"^) 1 + - + S ("*«)' + a ( (4,J "f + < 2 -» 

Here, i^ ln ^ and 3^ 1 ' 1 ^ are respectively the symmetric-trace- free (STF) mass and current multipole moments of the 
radiation field, and the antescript denotes the number of derivatives taken with respect to time. The luminosity is 
expressed as a function of time T, which represents proper time as measured by static observers at large distances. 
The total energy radiated, from initial time to final time Tf, is given by 

AE = / E{T)dT. (2.2) 

The radiative multipole moments in Eq. ( |2 . l| ) can be related to the multipole moments of the source. These are 
expanded in powers of v 2 ~ m/r to the order required for an expression for the luminosity accurate through 2PN 
order. They can be expressed as 
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The radiation measured at time T by an observer at a large distance R from the source is made up of two 
contributions. The first can be thought of as radiation propagating directly from the source; it depends on the state 
of the source at retarded time u ~ T — R (see below for a more precise expression). The second can be thought of 
as radiation scattered by the spacetime curvature as it propagates away from the source. This contribution, called 
the "tail," is emitted by the source at all times earlier than u. The tail term can also be decomposed into multipole 
moments. The dominant term, and the only one needed to evaluate E to 2PN order, is given by the following integral 
over retarded time along the source's trajectory C(it') prior to retarded time u Jl7[ : 
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(2.4) 



Here, u' parametrizes the trajectory up to the final point u' = 0. The constant iS is an arbitrary scale parameter 
which arises from the joining of two coordinate systems (see Sec. IIC for a discussion); we will show that our results 
do not depend on the value of this parameter. 

The last element needed to compute the energy flux are the equations of motion. For any two-body problem the 
relative acceleration a = a\ — a 2 in a harmonic coordinate system with origin at the center of mass, accurate to 2PN 
order, is given by jTsJ 
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We substitute Eq. ( |2 . 5 [ ) to the required order whenever time derivatives of the velocity appear in the various terms 
comprising Eq. (2.1). 
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Solving Eq. ( |2.5| ) yields the trajectory x(t). As in Newtonian mechanics, solving this equation is facilitated by 
making use of the first integrals: E the energy, and J the angular momentum. Both of these quantities are conserved 
to 2PN order. Below we will specialize to the case of a head-on collision, for which J = 0. We therefore only need an 
expression for the energy, which is given by 

f 1 o /"m\ 3, , . 1 , , /m\ 1 /m\ 1 /m\ 2 5 , „, E 

£ = /i{ 5 « 2 - ( 7 ) + g (1 " 3r?) , 4 + -(3 + „) ( 7 ) t; 2 + -77 (-) f 2 + - (-) + -(1 - 77? + 13 ^) v e 

+ 1(21 - 23 ?/ - 27 r, 2 ) (™) t; 4 + I77 (1 - 15 77) (™) i, 2 f 2 - jjr? (1 - 37?) (™) f 4 

+ 1(14 - 557/ + 4T7 2 ) + 1(4 + 69/?+ 12ry 2 ) (7) V - 1(2 + 15 tj) (^) 3 j. (2.7) 



B. Head-on collision 



In this subsection we use the equations listed above to calculate the gravitational-wave luminosity produced during 
the head-on infall and collision of two bodies of arbitrary masses. We consider two different situations. In the first 
(A), the infall proceeds from rest at infinite initial separation; in the second (B), the infall proceeds from rest at finite 
initial separation. 

There is only o ne d irecti on o f motion in this problem. We therefore have x = zn, v = in, r — z, and v = r = i. 
Using this, Eqs. (2.1) and (2.3) simplify considerably. Moreover, since the current moments are all proportional to 
J = 0, they all vanish. The mass multipole moments reduce to: 
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Here, the P s are STF products of a number q of unit vectors, as given by Pirani 120] and Thorne 16 1: 
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(2.8b) 
(2.8c) 
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(2.9b) 
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the P^'s are symmetric products of a number q of uni t ve ctors, and D symbolically represents the Kronecker delta. 
The compact notation to the right-hand-side of Eqs. (2.9), with which we represent fully symmetrized products of 
symmetric tensors by enclosing them within parentheses [e.g., (D ■ P^) or (D ■ D)], corresponds to expressions (2.1) 
and (2.2) in Ref. fll6fl - Notice that l^aii m ^+ (2.8a) also contains a factor P^; it arises via the term ( 2 )^pj m ^ ne 
integrand of Eq. ( |2.4[ ) . 

As the movement is directed along the z axis, the P^'s have components along z only (e.g., P^ — 5 lz 5 : > z ). The 
P s are therefore diagonal tensors, constant in time, and they factor out of the time derivatives in Eq. ([2.1|). As a 
result, we need only calculate the squares of the P s, 

2 p(3) p (3) _ 2 - (4) p (4) _ _8_ 
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and then insert them into the squares of ( n+1 )l !l in Eq. (2.1) 
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We must now evaluate the time derivatives of the ma ss m ultipole moments, substituting the expression for the 
acceleration z when needed. For a head-on collision, Eq. (2J5) reduces to 
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It is also helpful to make use of our expression (2.7) for the conserved energy in order to express z as a function of 
m/z. For a radial infall, we have 
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Defining zq to be the initial separation (at which z — 0), we can invert Eq. ( |2.12| ) for z(z, zo). For case (A), in which 
the initial separation is infinite, E(z) = E(oo) = 0, and we obtain 
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For case (B), in which the initial separation is finite, we have 
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We can then insert Eqs. ( 2.11| ), (2.13), and Eq. (2. IE) into the time derivatives of the mass moments, Eqs. (2.S), and 
calculate the gravitational- wave uminosity in terms of m/z and m/zQ. 

We discuss in the Appendix an interesting cancellation of terms which occurs in the case of an infall from infinity. 



C. Tail terms 



We now turn to the evaluation of the tail contribution to the gravitational- wave ene rgy flux. To this end, we must 
first consider the tail correction to the mass quadrupole moment, as given by Eq. (2.4). 

To leading order, we have l zz oc z 2 . Because z sw (m/z) 1 / 2 , each time derivative has the effect of multiplying 
this quantity by to 1 / 2 ;?" 3 / 2 . Integration with respect to retarded time then has the effect of multiplying the resulting 
expression by by z/z ~ z 3 / 2 to -1 / 2 . As a result, we have that -J oc (to/z) 3 ^ 2 \ zz , which shows that the tail correction 
is of 1.5PN order relative to the direct quadrupole term. Post-Newtonian corrections to Ixaii wm then be of 2.5PN 
order relative to the same term, and can be neglected in an expansion accurate through second post-Newtonian order. 
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the source's Newtonian trajectory. 
For case (A) — infall from infinity 



the tail term contributes to the luminosity via 
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where the relation between separation z and time t on the Newtonian trajectory is given by 
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The integral in Eq. (2.16) can be evaluated in closed form, and the result is 

71 5tt 



15 In 



6 V3 

For case (B) — infall from a finite distance zq — the tail term is 
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The Newtonian trajectory is here given by 
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where g(x) = y/xy/1 — x — arcsin \fx, with x = z/zq < 1. Evaluating the integral, we obtain 
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where Int(x) is given by 
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In this form, Int(z) has a very slow numerical convergence, because the argument of the logarithm vanishes at the 
lower limit. To evaluate it, we integrate by parts and change variables, so that 
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where y — arcsin y/x, and h(w) = sin w cos w — w. This expression is much better suited for numerical calculation. 

We now show that the gravitational-wave luminosity is independent of the arbitrary parameter S appearing in 
Eqs. ( |2.1g| ) and (2.21). For simplicity, we will only consider case (A) — infall from infinity; the argument for case (B) 
is entirely analogous. 

In the approach of Blanchet and Damour Jl7j ] , the energy flux at infinity is obtained by first solving the generation 
of gravitational waves in the near zone containing the source, using harmonic coordinates and expanding about flat 
spacetime, and then matching this solution outside the source to a solution of the vacuum equations in the far zone, 
expressed in radiative coordinates. In the near-zone solution, the waves propagate along "flat" null cones given by 



t — r = c, 

while in the far-zone solution, they propagate along the true null cones given by 

t — r — 2rnlnr = c , 



(2.24) 
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where c and d are constants. The two solutions are matched at an intermediate distance S ~ A ~ r/v, where A 
is the gravitational wavelength. This matching is the origin of the tail term, Eq. (2.18), together with its apparent 
5-dependence. 

When evaluating the energy flux, we must trace an observation made at time T by an observer at distance R back 
along the null cone to the corresponding time £cm (evaluated at r = 0, the center of mass) which determines the state 
of the source. The null cone that reaches the observer (T, R) is characterized by d = T — R — 2m\nR; the match at 
S gives d = c — 2m In 5, and the null cone at r = is characterized by icM = c. The result is 
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We find it convenient to rescale iS by m, and to define retarded time u = T — R — 2m\n{R/m). Then 
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This exhibits the explicit 5-dependence in the connection between z(t) and the observer. We now show that this 
5-dependence cancels the 5-dependence in the tail term. 
The trajectory at Newtonian order becomes 
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where z' = z[u); we have used the fact that tcM < 0, and have performed the expansions assuming m/|u| <C 1. The 
Newtonian contribution to ( 3 )T J is just the corresponding derivative of the quadrupole moment ij| = nz 2 P^ 2 \ 
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where the Newtonian expressions for i and 2 were substituted to yield a result in terms of z only. Applying the shift 
( [2.28 ) to this contribution, we obtain 
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The additional term in Eq. ( 2.30 ) cancels the dependence on ln(m/£>) in the tail term, Eq. ( 2.18 ). 

Thus the gravitational-wave luminosity is truly independent of the matching scale S. In practice therefore, we may 
assume that we have identified the orbital variable z at the correct retarded time, and eliminate the 5-dependence in 
the tail term by setting S = m. 



D. Luminosity and energy loss 



We now have all the elements need ed to calculate the grav itatio nal-w ave l uminosity for both types of collision, 
(A) and (B). We use the acceleration (2.11) and the velocities (2.13) and (2.15) to obtain the time derivatives of the 
multipole moments (2.8), including the tail terms, Eq. (2.18) and (2.21). We then square each term, inserting the 
squares of the STF tensors, Eq. (2.10). Finally, we compute E. 



1. Inf all from infinity 
For case (A) — infall from infinity — the luminosity is given by 
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In Eq. (2.31), within the curly brackets, we label the terms of order 0(1), 0(m/z), 0[(m/z) 3 ^ 2 ], and 0[(m/z) z ], as 
[N], [1PN]7[1.5PN], and [2PN], respectively. 
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In Fig. 1 we plot E/rj 2 as a function of z for different values of 77 = fx/m, from 77 — (test-body limit) to 77 = 0.25 
(equal-mass case). We see that E increases monotonically with increasing 77, and that the curves converge at large 
z; for such large separations, E/ij 2 ~ (16/15)(m/z) 5 , independent of r\. The plot also shows that E increases with 
decreasing z, until it reaches a maximum and then starts decreasing. This behavior signals the breakdown of the 
post-Newtonian expansion, as we now discuss. 



In Fig. 2a (for 77 = 0) and 2b (for 77 = 0.25) we separate the luminosity (2.31) into the various contributions 
corresponding to the [1PN], [1.5PN], and [2PN] terms, and plot their ratio to the leading-order, Newtonian [N] 
expression. Notice that both the [1PN] and the [2PN] contributions are negative, while the [1.5PN] contribution is 
positive, because of its logarithmic term. In both figures, for large z, the post-Newtonian terms are all smaller than 
the Newtonian expression. However, they are larger than might be expected a priori: at z = 100m, the [PN], [1.5PN], 
and [2PN] contributions would be expected to be of order 10~ 2 , 10~ 3 , and 10~ 4 , respectively. Instead, we find that 
they actually are of orde r 4 — 6 x 10~ 2 , 7 x 10~ 2 , and 4 x 10 -3 , respectively. This reflects the presence of large numerical 
coefficients in Eq. (2.31). As z decreases, the post-Newtonian terms increase until they become comparable to each 
other and to the Newtonian expression; this occurs within the interval 7 z/m < 10. For z <j 7m, the negative 
post-Newtonian corrections overcome the Newtonian term, and the luminosity formally changes sign. Evidently, the 
post-Newtonian approximation is no longer reliable in this region. Finally, comparing Fig. 2a with Fig. 2b, we see 
that the term which varies the most with 77 is the [1PN] contribution; the [2PN] term varies only weakly with 77, and 
the [1.5PN] term is independent of 77, as Eq. ( |2.3l| ) directly s how s. 

The total energy radiated, up to time 2/, is given by Eq. (|2.2|), 



AE(z f )= E(T)dT=- 



(2.32) 



The luminosity is given by (2.31) and z by Eq. (2.13). It is not possible to evaluate this integral in closed form. 
Instead, we approximate 1/z by its expansion in powers of m/z, up to the second order, and then integrate. We 
obtain 



AE(z f ) 




7 ( m 
5\/2 \Zf 



3/2 r 



53 5tt ir] m\ , 2 
h — = + 15 In — -hi - 
fi VS \ZfJ V3 



1 

66 



2 r 



81985 749 



3G 



-V 



18323 
16 



(2.33) 



In Fig. 3 we plot AE/r) 2 m as a function of Zf, for different values of 77. For large Zf, the energy radiated behaves as 
its Newtonian expression, AE/rfm ~ (16\/2/105)(to/z/) 7 / 2 , independent of 77. For small Zf, AE/rfm increases with 
increasing 77. In Fig. 3, the smallest value of Zf at which AE was calculated corresponds to where the gravitational- 
wave luminosity formally vanishes. 



2. Infall from a finite distance 
For case (B) — infall from a finite distance Zo — the luminosity is given by 



E = 
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(2.34) 



where x = z/zq and Int(x) is given by Eq. (2.22). This luminosity has characteristics similar to that calculated 
previously. The main difference is that here, contrary to case (A), the luminosity does not vanish at the beginning of 
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the infall. This, we shall now explain, is due to the fact that although the velocity vanishes at z = zq, the acceleration 
does not; that E does not vanish at z = zq merely reflects the time symmetry of the trajectory at t = t(zo). At 
Newtonian order the general expression for the energy flux, E = (8r; 2 /15) (m/r) 4 ( 12 v 2 — 11 f), necessarily vanishes 
at a moment of stationarity (v = r = 0). This is because it arises from an odd number (3) of time derivatives of 
a mass multipole moment. At higher order, however, contributions arising from even numbers of time derivatives 
appear, e.g., ( 4 )l yfc in Eq. (|2.l|) . These contain terms involving the acceleration, which does not vanish at zq. As a 
consequence, the luminosity also will not vanish. 

As in case (A), we see here also that the luminosity changes sign at some small value of z, when the negative higher- 
order terms become comparable to the positive Newtonian term. This signals the breakdown of the post-Newtonian 
method. 

The total energy radiated, from the moment of time symmetry to the final separation zt, is given by 

AE(z f ) = - H E(z)^, (2.35) 



where E(z) is given by Eq. ( 2.34| ) and z by Eq. (2.15). We have integrated this equation numerically for selected 
values of zo; the smallest value of Zf is chosen to be the one for which the luminosity formally vanishes. 

In Figs. 4a (for r\ = 0) and 4b (for -q = 0.25) we plot AE(zf)/i~i 2 m as a function of the final separation Zf, for 
different values of the initial separation zq, including the limit zo = 00 considered previously. We see that as zj 
decreases, the curves all approach each other. This results from the fact that most of the energy is generated near 
the end of the infall. 



III. HEAD-ON COLLISION OF A SMALL MASS AND A MASSIVE BLACK HOLE 

In this section we calculate the gravitational-wave luminosity produced during the radial infall, proceeding from 
rest at infinity, of a particle with small mass into a massive, nonrotating black hole. The smaller mass will now be 
denoted fi, and that of the black hole m. We assume fi <C to; in this limit fi and to are equivalent to the reduced 
mass and total mass of the previous section. 

This restriction on the mass ratio implies that the problem considered here is a limiting case of the one considered 
in the previous section, in which no restriction was put on the masses and the nature of the colliding objects. On the 
other hand, the problem considered here can also be seen to be an extension of the one considered in the previous 
section: contrary to Sec. II, we shall here put no restriction on the velocity of the infalling mass, and correspondingly, 
put no restriction on the strength of the gravitational field at the particle's location. While the results of the previous 
section were accurate through second post-Newtonian order, the results presented here are accurate to all orders in 
v/c. 

The stress-energy tensor associated with the infalling particle creates a small perturbation in the gravitational 
field of the nonrotating black hole (whose metric is given by the Schwarzschild solution) . Part of this perturbation 
represents the Coulomb field of the infalling particle; the remaining part represents radiative degrees of freedom, and 
these propagate away from the source as gravitational waves. 

We use the Teukolsky perturbation formalism |2l| to calculate E(T), the gravitational- wave luminosity as measured 
near future null infinity; the luminosity is expressed as a function of time T (to be defined precisely below). In Sec. IV 
this result will be compared to that obtained in Sec. II using post-Newtonian theory. 

The total amount AE — J E(T)dT of gravitational- wave energy produced during the radial infall of a particle into a 
Schwarzschild black hole was first calculated by Davis, Ruffini, Press, and Price H. This calculation was subsequently 
generalized to nonradial motion by Detweiler and Szedenits Q . Infall into a Kerr black hole was considered by Sasaki 
and Nakamura Q (radial motion) and Kojima and Nakamura (nonradial motion). 

The Teukolsky perturbation formalism is described in detail in Ref. we shall make frequent use of equations 
contained in this paper. In the Teukolsky formalism, the gravitational perturbations are described by a single, 
complex-valued function ^4, which represents a particular component of the perturbed Weyl tensor. The second- 
order differential equation governing 'J 4 can be completely separated by decomposing ^4 into Fourier modes and 
spherical-harmonic components. The radial function Eim(w;''), where lo is the angular frequency and L, M the 
spherical-harmonic indices, satisfies an inhomogeneous, second-order, ordinary differential equation, whose source 
term is constructed from the particle's stress-energy tensor. This equation can readily be solved by means of a 
Green's function. 

As shown in Sec. II of Ref. |2^] , the gravitational waveforms, h + and h x , as measured by a fictitious detector situated 
near future null infinity, can be obtained from the asymptotic behavior Rlm{w, r — > 00) of the radial function. With 
a slight change in notation, Eq. (2.14) of Ref. |2^] becomes 
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h(T,R,e,<t>) = ^$^m(«)- 2 }lm(M)- (3.1) 

LM 

Here, h = h + — ih x , T is proper time at the detector's location, R is the distance from the source to the gravitational- 
wave detector, and 

u = T-R* =T - R-2m\n(^j- -l^j (3.2) 

is retarded time; T and R are Schwarzschild coordinates, which are distinct from the harmonic coordinates used in 
the preceding section. The functions Zlm(u) represent the multipolc moments of the gravitational- wave field; they 
are expressed as the Fourier integrals 

Zlm{u) = J Z LM (Lu)e^ u duj, (3.3) 

where Zlm(u) will be given below. Finally, -2Ylai(0, (f>) are the spherical harmonics of spin- weight —2 |g3|] . 
In terms of h given above, the gravitational-wave luminosity can be expressed as 

R^CdhOh 

v ' 16tt / dTdT ' y ' 



where an overbar denotes complex conjugation, and the integration is over solid angles. Substituting Eq. fl3.1| ) into 
( |3.4[ ), and using the orthonormality of the spin- weighted spherical harmonics, we obtain 

where 

Zlm(u)= I (-iuj)Z LM (uj)e-^ u du;. (3.6) 



477 

LM 



It is straightforward to follow the prescription described in Ref. |2^] and to calculate an expression for Zlm(u) for 
the problem under consideration. The first step is to describe the motion of the particle, whose world line is taken 
to be a marginally bound, radial geodesic of the Schwarzschild spacetime. We take E = u t = 1, where u a is the 
particle's four-velocity, and 9 — along the world line. The geodesic equations are 

dr /2rn\ 1/2 dt 1 



dr \ r J ' dr / ' ^ 

where r denotes proper time and / = 1 — 2m/ r. These imply the following relation along the world line, where 
x = (r/2m) 1 / 2 : 

t(r) = -Amg(x), (3.8a) 

q(x) = -x 3 + x + - In — -. (3.8b) 
yw 3 2 x + 1 y ' 

Proceeding along the lines described in Sec. II of Ref. ^2|, we obtain the following expression for Zlm(ui): 

Zlm(u) = Z l (uj) S M o- (3.9) 

The fact that only the modes with M = contribute to the radiation reflects the axial symmetry of the problem. We 
also obtain 



2muj (L - l)L(L + 1)(L + 2) - Ylimu V ^ A%{u) J 2m (i + ^^fr^ 

(3.10) 
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The quantities A™(uo) : Tl(ui), and Xl(lo\t) have not yet been introduced. We shall explain their meaning in the 
following paragraph. 

The function Xl(uj; r) is a solution to the Regge- Wheeler equation p4|, 



dr* 2 



+ u 2 + f 



L(L + 1) 



6m 



>X L (u-r) = 0, 



(3.11) 



where d/dr* = fd/dr. It is chosen so as to have the asymptotic behavior 

X L (uj;r -> 2m) - e~ Mr * 
near the black-hole horizon. Correspondingly, 

X l (cj- r -f oo) ~ ^» e — * + 0(e^*); 



(3.12) 



(3.13) 



this equation defines the constant AjP(cj) appearing in Eq. (3.10). Our expression for Zl(u>) also involves the first-order 
differential operator 



T l (lj) = 2(1 - 3m/r + iur)rf — + f [L(L + 1) - 6m/r] 

ar 

+ 2iwr(l — 3m/r + iwr). 



(3.14) 



The strategy to calcul ate th e gravitational- wave luminosity is the following. For given val ues o f L and u>, the 
Regge- Wheeler equation ( 3.11 ) is integrated numerically, starting near r — 2m and using Eq. ( 3.12 ) to specify the 
initial conditions. The Regge- Wheeler function and its first d erivative are evaluated at values of r lyin g in the interval 



(2m, oo). The constant A™(oj) is computed using Eq. j3.13j ). The integral to the right of Eq. ( 3.1C ) is then carried 



out, numerically, for the given values of L and w. These steps are repeated for many relevant values of the fre quen cy, 
and the integral to the right of Eq. ( |3.6| ) is evaluated for many values of u. Finally, the sum over L in Eq. (3.5) is 
carried out, and E(u) is obtained for t he selected values of u. 

The integral to the right of Eq. ( 3.1 0| ) is actually divergent, and must be regularized before carrying out the program 
outlined in the preceding paragraph. To do this we follow Detweiler and Szedenits 0], who showed that the divergent 
behavior can be removed by performing an integration by parts, and then discarding a (formally infinite) boundary 
term. The ultimate justification for this somewhat dangerous procedure comes from the eventual agreement with the 
previous analysis of Davis et al. pi , which is based on a manifestly finite perturbation formalism. 

We denote the integral to the right of Eq. ( 3.1 0| ) by Il(u>). It is easy to see that this integral can be split into 



convergent and divergent parts, so that I = I c , 



Idiv 



/div (indices are now suppressed for simplicity), with 
a{r) e tuJt( - r) £X(r)dr, 



where C = fd/dr + iu>, and 



a{r) 



2r v /r/2r 



(072771 + 1) 



(1 — 3m/r + iur) ; 



(3.15) 



(3.16) 



the remainder of the integral gives J conv . 
To regularize /div we rewrite Eq. ( |3.15 ) as 



^div — 



ae luJt £X + —(be wt CX 
dr V 



2m 

be iuJt CX 



dr 



2m 



(3.17) 



and seek a function b(r) such that the new integral converges. The boundary term at infinity will then be seen to be 
infinite, and will be discarded; the boundary term at the horizon will be seen to vanish. It is easy to check that if 
b(r) satisfies the differential equation 



db f dt 1 j 

dr \dr f , 



a = 0, 



(3.18) 
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then the new integral will converge. A particular s olutio n to Eq. (3.18) is b(r) = — 4imo-'~ 1 (l + :r) _1 (l + x + 2imujx 3 ) 
where x = (r /2m) 1 1 2 . Substituting this int o Eq. ( [3.17 ), discarding the boundary terms, and combining with J conv 
we find that the integral to the right of Eq. ( 3.1 0| ) can be written in the regularized form 



1 + imuj (r/2m) 3 / 2 



2m imuj(r/2m) 2 
x [L(L + 1) - 6m/r] X l (uj; r) dr. 



(3.19) 



It is evident that this integral is convergent at r — oo. 

The numerical code written for the purpose of calculating E(u) was built upon FORTRAN subroutines given in 
Numerical Recipes |2E) |. The Regge- Wheeler equation is integrated to high accuracy using the Bulirsh-Stoer method, 
which is also used to evaluate the integral ( B.19| ). [This integral is truncated at a radius r large compared with 
1/lu; the remainder is calculated analytically using the asymptotic form for the Regge- Wheeler function, as given in 
Eq. ( 3. 13 ).] These steps are repeated for many values of oj; w e have calculated more than 600 points in the relevant 
interval < mui < 0.9 The Fourier transform of Eq. ( |3.6[ ) is carried out by first selecting a value for u, and then 
evaluating the integral using a Romberg integrator. An interpolator is used to compute the integrand at those values 
of to which the integrator selects; these are in general distinct from the values selected when previously calculating 
Zl(ui). These steps are repeated for many values of u, lying in the relevant interval —80m < u < 40m. Finally, the 
sum over L is carried out. Due to the rapid convergence of the multipole expansion, only the modes with L — 2 and 
L = 3 need be included in the sum. 

A plot of E(u), the gravitational-wave luminosity as a function of retarded time, is displayed in Fig. 5. It shows 
two main features: bremsstrahlung radiation and black-hole quasi-normal ringing. The bremsstrahlung radiation 
dominates the luminosity at retarded times u < —10m. This radiation is generated by the particle's infalling motion at 
early times (large distances from the black hole) , and propagates directly from the particle to the detector (scattering 
by the spacetime curvature is insignificant). The black- hole ringing dominates the luminosity at retarded times 
U > —10m. This radiation is generated by the dynamics of the gravitational field in the vicinity of the black- hole 
horizon; it is scattered many times by the spacetime curvature, and is essentially disconnected from the particle's 
motion. It is clear that quasi-normal ringing dominates (by approximately two orders of magnitude) the energetics 
of the problem. 

By integrating E(u) over retarded time we were able to reproduce the standard result for the total energy radiated, 



AE= (0.0092 + 0.0011 



)^ 2 /m, 



(3.20) 



first calculated by Davis et al. [Q. In Eq. ( 3. 20] ), the first term comes from the L — 2 mode, while the second 
corresponds to L = 3. The total coefficient, as computed in Ref. M and reproduced here, is equal to 0.0104. 



IV. COMPARISON BETWEEN POST-NEWTONIAN AND PERTURBATION-THEORY RESULTS 



The purpose of this section is to compare the results obtained in Sec. II using post-Newtonian theory to those 
obtained in Sec. Il l usin g black-hole perturbation theory. More precisely, we shall consider the r\ — > limit (where 
r\ = /i/m) of Eq. ( 2.31 ) and compare it to the numerical results displayed in Fig. 5. While the numerical results 
are valid to all orders in the post-Newtonian expansion, the analytic expression is accurate only through second 
post-Newtonian order. 



In the limit of small mass ratios, Eq. (2.31) takes the form 



^,/„s 16 ,/ m\ 43 ( m \ — 



71 57r / m 

TT + 73 + 151n U 



Mu,§ 
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(4.1) 



For convenience, we have made slight changes to the notation used in Sec. II: we have replaced z by r^, where the 
subscript indicates that the coordinates are harmonic. The gravitational-wave luminosity is expressed as a function 
of time T/j, which operationally has the same meaning as in Sec. Ill: it corresponds to proper time as measured by a 
static gravitational-wave detector situated near future null infinity. 

Equation (4.1) is incomplete without the relationship between Th (observer harmonic time at large distances), 
Which appears to the left, and r/j, which appears to the right. This relationship is discussed in Sec. IIC: If rh{th) 
describes the particle's world line in harmonic coordinates, then the quantity Th appearing to the right of Eq. fl4.1| ) is 
rh(u), where u is the retarded time 
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Th — Rh 



o 1 ( Rh 
2m In — 

V TO 



(4.2) 



Here, Rh ^ to is the distance from the detector to the system's center of mass (wh ich here is identical to the fixed 
position of the large mass to). Notice that Eq. ( [1.2] ) ta kes the same form as in Eq. ( |3.2| ), in the limit of large R/m, 
apart from two slight differences. The first is that Eq. (3.2) involves Schwarzschild coordinates, T$ and Rs, instead 
of harmonic coordinates. The relationship between these coordinates is 



const. 



Rs = Rh 



(4.3) 



The second is that in Eq. ( |4.2| ) Rh is divided by to, while in Eq. ( |3.2| ) -Rs is divided by 2m. This difference originates in 
the freedom in selecting the value of the parameter S, as was discussed in Sec. IIC. This corresponds to the freedom 
in choosing the origin of the retarded time u. In particular, the origin of time might be chosen differently in the 
post-Newtonian and perturbation-theory calculations. 

Because of this ambiguity in the origin of u, it would be unwise to compare directly Eq. (4.1) with Fig. 5. A much 



better way of carrying out the comparison is to re-express the luminosity function in terms of a fully unambiguous 
parameter. We shall choose to parametrize the world line with the Schwarzschild coordinate rs, and express the 
luminosity in terms of this parameter. It might appear that this re-expression has already been effected (apart from 
the transformation rs = rh + to) in Eq. (3.2). We shall see, however, that the right-hand side of Eq. ( |3.2[ ) is not the 
desired expression. 

Let us begin with the perturbation-theory results for E(u), and let us adopt a mapping between u, the time 
coordinate running along future null infinity, and rs, the adopted parameter along the particle's world line. The 
simplest and most natural mapping is depicted in Fig. 6: we connect the event labeled by rs on the world line and the 
event labeled by u on future null infinity with a radial, outgoing null geodesic. This null geodesic is taken to propagate 
in the direction directly opposite to the direction in which the particle moves. Radial, outgoing null geodesies are such 
that ts — rs — 2mhx(rs/2m— 1) is constant along th em. The relation u(rs) therefore follows directly from Eq. fl3.2|), 
with rs substituted for R, and ts(rs), given by Eq. (3.8), substituted for T. The result is 



u(rs) = —4m 



^x 3 + \x 2 + x + ]n(x - 1) 



(4.4) 



where x = (rs /2m) 1 / 2 . Using Eq. (4.4), the luminosity function can be plotted as a function of rs, the parameter 
along the world line. This plot is displayed as a solid curve in Fig. 7. 

The figure shows that the numerical curve tends to become pointwise unreliable as the radius increases beyond 
rs = 20to: the computed curve undergoes small oscillations about a mean curve which presumably represents the 
true luminosity function. These oscillations are due to numerical error. More precisely, they are a consequence of the 



fact that the integration over angular frequencies, cf. Eq. (3.6), must necessarily be cut off at some finite upper bound 
<x>o. Reducing this upper bound produces larger oscillations about approximately the same mean curve. 

The mapping between world line and future null infinit y ad opted in Eq. (4.1) is different from the one constructed 
in the preceding paragraph. This is understood from Eq. (4.2), which involves the distance Rh to the system's center 
of m ass, instead of the distance to the particle itself. Therefore our plot for E(rs) cannot be compared directly with 
Eq. (4.1). To carry out a meaningful comparison, we must first map rA, the position of the particle at the retarded 
time of the center of mass, into the true-retarded position rs, as illustrated in Fig. 8. (The coordinates used here are 
the harmonic coordinates, but we suppress our use of the subscript h.) We now turn to this task. 

Figure 8 illustrates how to construct the mapping from rA to r B . As mentioned previously, the coordinates used 
for this purpose are the harmonic coordinates. Let t = p(r) represent the particle's world line, and let t = q(r) + u 
represent a light ray propagating in the opposite direction and meeting future null infinity at retarded time u. The 
center of mass is assumed to have the fixed position r = 0. The radius rA is determined by inverting the equation 

p(r A ) =q(0)+u. (4.5) 

The radius rs, on the other hand, is obtained by solving 

p{r B ) = q(r B ) + u. (4.6) 

By eliminating u from these equations, we obtain the desired (implicit) relationship between rA and r B : 

p{r A )=p{r B )-q{r B ) + q{0). (4.7) 



Our goal now is to transform Eq. (4.7) into something more concrete 
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We need to derive appropriat e exp re ssion s for p(r) and q(r). The situation for q(r) is quite simple. An expression 
for it can be derived from Eqs. ( 2.24| )-( 2.27 ). With S consistently set equal to m, we obtain simply 



q(r) = r. 



(4.8) 



A n exp ression for p(r) can be derived by integrating dr/dt, which was written down in Sec. IIB. In the limit 77 — s> 0, 
Eq. (2.13) becomes 



dr 
~dt 



in- \ii I * 5 27 o 

— = -y 1 1 — y h — y 



(4.9) 



where y — 2m/r; Eq. (4.9) is valid through second post-Newtonian order. Integrating Eq. (4.9) would give an 
expression for p(r) also accurate to second post-Newtonian order. This would imply that the mapping is carried out 
with the same degree of accuracy as the calculation of E, as giv en in Eq. (4.1). However, we have here the opportunity 



of using an expression for dr/dt that is more accurate than Eq. (4.9). This comes about because in the 77 — > limit, the 
equations of motion are given exactly by the Schwarzschild expressions (3.7). Once written in terms of the harmonic 
coordinates, these equations imply 



dr 
~dt 



fL = _ w l/2 



-3/2 



(4.10) 



Expanding this in powers of y pq| , we obtain 



dr 
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(4.11) 



which generalizes Eq. (4.9). For the purpose of carrying out the mapping we will deal with Eq. (4.11), whose integration 
gives p(r): 
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The mapping between r^ and rg is given by Eqs. (4.7) and ( |4.8| ): p{rA) — p{^b) ~ r B- This equation could be 
inverted numerically to yield r^rs). However, the following approximate, analytic inversion proves t o be sufficiently 
accurate for our purposes. First, we define an auxiliary variable Z = [— Zp{r a) / ~ 2 / 3 . Next, Eq. ( 4.12Q is used to 
express Z as a power series in yA = 2m/rA\ this series is then inverted to yield yA as a power series in Z: 
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Finally, we use our mapping equation to write Z 
Z in terms of ys = 2m/rB'- 



-3p(r b) I '4m 4- 3r b I '4m] 2 ^ 3 ; substitution of Eq. (4.12) then gives 



Z = y B [l 
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(4.14) 



We have obtained the desir ed inversion: a value of rs is selected, Z is calculated usin g Eq. ( 4.14 ), and is 
then obtained from Eq. ( 1.13j ). It is this radius which must be substituted into Eq. (4.1) in order to obtain the 
gravitational- wave luminosity. This can then be plotted against the Schwarzschild radius r$ = rs + m, which gives 
the dotted curve in Fig. 7. 

The accuracy of the post-Newtonian expression can be ascertained from the figure. We find that the post-Newtonian 
value differs from the numerical one by a factor of approximately 2.5 at r$ = 10m, 1.5 at r$ = 20m, and 1.3 at 
rs = 30m. We see that the post-Newtonian values tend to converge at large distances to the exact values, but that 
the rate of convergence is slow. We will explore this issue of the convergence of the post-Newtonian values more fully 
in the following section. 
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V. ACCURACY OF THE POST-NEWTONIAN EXPANSION 



We have se en in the preceding section that the post-Newtonian expression for the gravitational-wave luminosity, as 
given by Eq. ( |4.lD , only approximately reproduces the exact, numerical curve, as represented in Fig. 7. And indeed, 
the degree of accuracy is far less than could be expected: An expression valid to second post-Newtonian order could 
be expected to have a fractional accuracy of order (m/r) 3 ; for r = 10m this is of order 10~ 3 . Instead, comparison 
with the numerical results shows that the fractional accuracy is actually of order 0.5. 

In the previous section, the comparison between the post-Newtonian and perturbation-theory results could be 
carried out for r$ < 30m only. The purpose of this section is to push the comparison to larger values of the radius, 
in spite of the fact that the numerical results are not reliable beyond r$ — 30m. 

To this end, we consider the following calculation, based on the hybrid formalism proposed by Kidder, Will, and 
Wiseman [[l9| . We calculate the gravitational- wave luminosity associated with the radial infall of a particle with small 
mass into a much more massive, spherically symmetric object (not necessarily a black hole). The masses are denoted fi 
and m respectively, and we assume /i <C m. We do so by usi ng th e multipole expansion of Eq. (2.1), and by employing 
the exact, Schwarzschild equations of motion given by Eq. ( |4.10| ). Taking advantage of the rapid convergence of the 
multipole expansion, only the mass quadrupole term will be considered; tail effects will also be ignored. 

The calculation proceeds as follows. In a Cartesian coordinate system based on the harmonic coordinates, the 



system's trace-free quadrupole moment can be expressed as V 3 
Taking three time derivatives, we obtain 



fir 



,(2) 



where P 



(2) 



<n\ / d 3 r dr d 2 r\ 

()iJ =H r ^ +3 ^> 



,(2) 



diag(-l/3,-l/3,2/3). 



(5.1) 



Use of Eq. ( 4.1 0| ) and substitution into the first term of Eq. (2.1) finally yields 

16 



E =15^ 



TV 



(5.2) 



where 
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(5.3) 



Equation (5J2) with TZ — 1 gives the standard quadrupole- formula result; the factor TZ 2 is the relativistic correction. 
Notice that TZ te nds toward unity at large r, and that TZ = at r = m (the e vent horizon) and r = 3m. 

Equation ( |5.2| ) is to be interpreted as the fully relativistic analogue of Eq. (LI). This expression can be compared 
with the exact, numerical curve of Fig. 7, provided that the mapping illu strated in Fig. 8 is carefully carried out. 
The result is displayed as the dotted curve in Fig. 9. We see that Eq. (5.2) reproduces the numerical results to very 
high accuracy for r$ > 7m. In particular, this expression does much better than its post-Newtonian analogue. This 
agreement sho uld not be considered to be of deep physical significance; it is essentially a fluke. However, it will allow 
us to use Eq. (5.2) as a tool for determining the large-distance accuracy of the post-Newtonian results, as we now 
explain. 

Although the numerical curve cannot, because of numerical error, b e re liably constructed for radii larger than 
approximately 30m, it is most plausible that the agreement with Eq. (^^) would persist at larger values of r. It 
therefore appears appropriate to investigate the large-distance convergence property of the post-Newtonian luminosity, 
with respect to the exact curve, by comparing Eq. (4.1) directly with Eq. ( |5.2[ ). This comparison, which has the major 
advantage of not involving the mapping illustrated in Fig. 8, reveals that the relative difference between the post- 
Newtonian values and the exact ones is 20% at ~ 45m, 10% at ~ 75m, and 1% at ~ 400m. These results 
establish the slow convergence of the post-Newtonian expansion. 



ACKNOWLEDGMENTS 



Part of this work was carried out with the aid of the symbolic manipulator MACSYMA. This research was supported 
in part by the National Science Foundation under Grant No. PHY 92-22902 and the National Aeronautics and Space 
Administration under Grand No. NAGW 3874. 



15 



APPENDIX: CANCELLATION OF POST-NEWTONIAN TERMS 



The purpose of this Appendix is to point out a remarkable cancellation of terms occurring in the calculation of the 
gravitational-wave luminosity in the case of a head-on collision proceeding from rest at infinity. 

Substituting Eq. ( 2.13] ) into Eqs. (2.8) we find that the mass multipole moments can be expressed as expansions in 
powers of m/z only. To second post-Newtonian order, the generic form is 



^ zz 

TZZZ 
tZZZZ 



\XZ 



An 



A 



3,0 



■A 2A 
A 3<1 
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(Ala) 
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= fj,z 4 A ij0 . 



(Ale) 



Assuming that such an expansion holds at least to order (m/z)™, we write the mass n-pole moment in the form 



=^z n {A nfi + A n , x [- 



Ar. 



(?)' 



O 



m\ n + 1 

z 



(A2) 



We have represented the first few moments in Fig. 10, in which a row corresponds to a given multipole order, and a 
column to a given power of m/z in the expansion. The thick dashed line separates the coefficients needed in order to 
express the luminosity accurately through 2PN order; these appear to the left of the boundary. 

First, we observe that in Eq. (|A2|), the term of order 0[(m/z) n ) is just the constant A n n [im n . The contribution 
from this term to the energy flux therefore vanishes when time derivatives are taken. We designate this cancellation 
in Fig. 10 by the darkly shaded diagonal beginning at ^2,2- 



Second, because z oc z 1 / 2 to leading order, the k th time derivative of a term of the form Ai 



with I > n, 



goes like z ( l - n - 3k / 2 ) , Consequently, if I — n = mo d 3, this term becomes a constant after k = 2(1 — n)/3 derivatives, 
and will vanish after an additional derivative is taken. But since each Z-pole moment is differentiated I + 1 times to 
calculate the energy flux, this term will vanish whenever fco < I — 1, or I + 2n > 3. Such terms are indicated in Fig. 10 
by the lightly shaded diagonals. 

This does not imply, however, that contributions from the PN order corresponding to a shaded region do not appear 
at all in the time-derivatives of the corresponding STF moment. If the term in question is a higher-order correction to 
the moment (e.g., ^2,2 or A^i), then contributions to ( 3 )i lJ of the same order could be generated by correction terms 
in the equations of motion, applied to derivatives of lower-order terms. Only if the term is a leading-order term (e.g., 
A^fi or ^6,0) is the contribution at that order identically zero. It is interesting to note that because of this effect, the 
coefficient A2.2 is not explicitly needed in calculating the energy flux to 2PN order. This coefficient has only recently 
been calculated pj| . 

In the case of infall from a finite distance, the argument used above fails. To leading order, the velocity is now 



m 



1/2 



(A3) 



and the 2PN expansion for i lJ contains powers of both (m/z) and (m/ zq). For example, terms of the form m 2 / zzq 
and (m/zo) appear in the 2PN expansion for the quadrupole moment; when multiplied by fiz 2 , these terms are not 
killed by time derivatives. 
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FIG. 1. Gravitational-wave luminosity, in units of rf , as a function of separation z, for case (A): infall from infinite initial 
separation. We show the luminosity curves for n = {0, 0.15, 0.25}. 

FIG. 2. A plot of the absolute value of the 1PN, 1.5PN, and 2PN contributions to the luminosity divided by the Newtonian 
expression. The 1PN and 2PN contributions are negative; the 1.5PN contribution is positive. Fig. 2a shows the test-body limit 
(n = 0), and Fig. 2b shows the equal-mass case (rj = 0.25). 

FIG. 3. Total energy radiated during an infall from infinite initial separation to final separation Zf. We show curves for 
rj = {0,0.15,0.25}. 

FIG. 4. Total energy radiated during an infall from initial separation zo to final separation Zf, including the limit zo = co. 
We show curves for zo/m — {15, 30, 50, oo}. Fig. 4a shows the test-body limit (r/ = 0); Fig. 4b shows the equal-mass case 
(r) = 0.25). 

FIG. 5. A plot of E /rf (on a logarithmic scale) as a function of retarded time u/m, as calculated using perturbation theory. 
The luminosity represents bremsstrahlung radiation for u/m < —10 and black-hole quasi-normal ringing for u/m > —10. 
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FIG. 6. Conformal diagram representing the Schwarzschild spacetime. Shown are: past future null infinity (lower diagonal 
with positive slope), future null infinity (higher diagonal with negative slope), the future horizon (higher diagonal with positive 
slope), the past horizon (lower diagonal with negative slope), and the singularity (broken horizontal). We also represent the 
world line of the infalling particle, parametrized by the Schwarzschild coordinate rs, and the radial, outgoing null geodesic 
reaching future null infinity at retarded time u. 

FIG. 7. A plot of E/rj 2 (on a logarithmic scale) as a function of rs/m, the parameter along the particle's world line. 
The solid curve represents the exact, numerical results obtained using perturbation theory. The dotted curve represents the 
post-Newtonian approximation given by Eq. (4.1). 

FIG. 8. The mapping between va, the center-of-mass-retarded position, and rs, the true-retarded position. The curve 
t — p(r) represents the world line of the particle. The curve t = q(r) + u represents an outgoing null ray propagating in the 
opposite direction, reaching future null infinity at retarded time u. The center of mass is located at the center r = 0. 

FIG. 9. A plot of E/rj 2 (on a logarithmic scale) as a function of rs/m, the parameter along the particle's world line. 
The solid curve represents the exact, numerical results obtained using perturbation theory. The dotted curve represents the 
quadrupole-relativistic approximation given by Eq. (5.2). 

FIG. 10. Cancellations of terms in the case of infall from infinity. Each row represents a multipole order, with a common 
factor z % , while each column labels the contribution of order 0(m/z) J to the multipole moment in a post-Newtonian expansion. 
Tail terms are ignored. A consistent expansion of the luminosity up to a certain order n involves the coefficients included in 
the upper-left corner of the table (the requirement for an expansion of the luminosity accurate to 2PN order is marked by the 
heavy dashed line) . The coefficients which have currently been calculated explicitly are indicated in bold face. 
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